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Abstract 

Testing the stability of the electroweak vacuum in any extension of the Standard Model 
Higgs sector is of great importance to verify the consistency of the theory. Multi-scalar 
extensions as the Minimal Supersymmetric Standard Model genetically lead to unstable 
configurations in certain regions of parameter space. An exact minimization of the scalar 
potential is rather an impossible analytic task. To give handy anal}Tic constraints, a spe¬ 
cific direction in field space has to be considered which is a simplification that tends to 
miss excluded regions, however good to quickly check parameter points. We describe a 
yet undescribed class of charge and color breaking minima as they appear in the Mini¬ 
mal Supers}nnmetric Standard Model, exemplarily for the case of non-vanishing bottom 
squark vacuum expectation values constraining the combination in a non-trivial 
way. Contrary to famous A-parameter bounds, we relate the bottom Yukawa coupling 
with the supersymmetry breaking masses. Another bound can be found relating soft 
breaking masses and p only. The exclusions follow from the tree-level minimization 
and can change dramatically using the one-loop potential. Estimates of the lifetime of 
unstable configurations show that they are either extremely short- or long-lived. 
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1 Introduction 


A complete analysis of the vacuum structure in any quantum field theory needs a considera¬ 
tion of the effective potential to all orders which is more than an honorable task. Important 
contributions to the effective potential in the Standard Model and supersymmetrized ver¬ 
sions at one and way more loops have been (partially) determined [1-6]. The more loops 
the more difficult is also the task to find the global minimum which shall determine the 
vacuum state of the theory. Numerical solutions to that problem exist in the Minimal Su¬ 
persymmetric Standard Model (MSSM) where both the effective potential as well as the 
(expected-to-be) global minimum are calculated and determined purely numerically [7,8]. 
Supersymmetry (SUSY) generically tends to stabilize the potential as negative fermionic 
loop contributions are compensated by the corresponding bosonic ones. The superpartner 
spectrum on the other hand brings additional directions in scalar field space that potentially 
invalidate the electroweak Higgs vacuum at the classical level. A physical viable supersym¬ 
metric extension has to take care of the additional parameters in a way that the “desired” 
vacuum is the true vacuum of the theory. 

The consideration of the one-loop effective potential, which can be very efficiently done 
via the famous formula of Coleman and Weinberg [1], leads to a first understanding of 
non-trivial minima. We have 
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where the sum runs over all fields f in the loop and counts gauge degrees of freedom 
like Cquark = 3 (spiu degrees of freedom are covered by the supertrace STr). The field- 
dependent mass eigenvalues A4^(0) are generically the eigenvalues of the Hessian matrix 
of the full scalar potential and the field (p represents any t37pe of scalar field value which 
is still present in the masses (do not set remnant field values to zero, they correspond to 
vacuum expectation values (vevs) at local or global minima of the potential). Additionally, 
there is a polynomial Pf^cp) which is renormalization scheme dependent and in the most 
common cases a constant. The renormalization scale is given by Q. 

The one-loop potential is known to develop an imaginary part [8-11] which is of no im¬ 
portance in the discussion of tunneling times from false to true vacua but opens the access 
to non-standard vacua; an imaginary part in the one-loop effective potential is related to a 
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non-convex tree-level potential at that point. ^ A non-convex potential means that the sec¬ 
ond derivative is negative which corresponds to a tachyonic mass eigenvalue < 0. 

The tachyonic mass, however, would only be present at the minimum (which by definition 
is locally convex). So, the existence of a non-convex direction points towards a minimum in 
that direction unless the potential is unbounded from below, which would be even worse. 
Finding the critical field value at which the non-convex direction opens is trivial as we shall 
see. The question is rather whether the non-standard minimum is deeper than the stan¬ 
dard one and therefore allows for a vacuum-to-vacuum transition which can be figured out 
analytically under certain circumstances. 

We first consider the loop corrected Higgs potential in the MSSM including SUSY loop 
contributions from the third generation (s)fermions. The tree-level part is given by the mass 
terms and the self-couplings which are gauge couplings. The one-loop part is given by the 
logarithms of Eq. (1) which also follow from the direct calculation [11]. We borrow the 
notation from [11] and define the effective potential as 


VeS = Vo + v! + V; + vl’ + V^^ 
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where the abbreviations x^ ^ and y^ ^ are 
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The soft SUSY breaking masses enter as m^, and and we defined Ml^ = + 

mlh)/2. The trilinear soft breaking couplings in the up and down sector are given byA^ 
and A^,, respectively. Yukawa couplings are denoted as 7^^, and jU is the jU parameter of 
the superpotential in the MSSM. The mass parameters of the tree-level Higgs potential are 

^It is actually related to a branch point of the logarithm in Eq. (1) that appears for a zero mass eigenvalue. 
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+ ImP, rnl' 2 ^ = + li^l^ 01 ^ 2 '^'" = with the soft breaking masses and 

for the and doublet, respectively; is the soft breaking bilinear term 
We consider only third generation superfields which couple with large Yukawa couplings to 
the Higgs doublets: 

W = /iHd-H, + Y,H,-QLfR-YbHd-QA- (4) 

The left-handed doublet field is Ql = (TljBl) and the two Higgs doublets = (h+, h°) and 
Hd = (h^, —hd); SU(2)L-invariant multiplication is denoted by the dot-product. The SU(2 )l 
singlets are put into the left-chiral supermultiplets = {t*, t^} and B^ = {b*, b^} with the 
charge conjugated Weyl spinors and b^. 

The effective potential of Eq. (2) obviously develops an imaginary part beyond the branch 
point of the logarithms ln(l ±x -l-y). We want to give a physical meaning of this branch point 
without reference to an imaginary part of the effective potential, since ^ In ((l±x-l-y)^) does 
not reveal any imaginary part—nevertheless, this logarithm gets singular where ~ J = 1 
though the potential itself stays finite. This point determines (for fixed parameters) a critical 
Higgs field value for which one mass eigenvalue gets tachyonic. The effective potential is 
a function of the (classical) field values which correspond to vacuum expectation values at 
the minimum. In the direction of the negative mass square, the potential drops down and 
therefore develops a CCB vacuum. 

Moreover, for certain parameters, the potential of Eq. (2) develops a second minimum 
in the direction of a standard Higgs vev which always lies beyond the branch point of one 
of the logarithms [11]. Expanding around this second minimum, one finds exactly one 
negative sbottom mass square (in the region of large jU and tan^) which hints towards a 
global minimum including a sbottom vev. The second minimum as depicted in [11] is an 
artifact of holding bL R = 0: the global minimum lies at a point with both (bL ^) 7 ^ 0 and 
K) 7 ^ V,. ’ 

We take the existence of the critical field value serious and first figure out its meaning 
for the development of such a CCB minimum. For simplification we now restrict ourselves 
in the following to (t^) = (t^) = 0 and also do not consider stau vevs. Let us consider for the 
moment a fixed value of the down-t 3 q 3 e Higgs field, b° = v^ and set A^, = 0. The critical field 
value is then obtained by solving — y^, = 1 with x^, and y^, given in Eq. (3) : 
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with and jU, as well as the Higgs field assumed to be real. The bottom 
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Yukawa coupling suffers from SUSY threshold corrections and reads Y^, = 
with including the Higgsino corrections ~ jjA^ tan P [12-15], which can be dominant 
over the gluino-induced threshold correction for large /r tan ^ and large gluino mass. Both 
gluino and higgsino contributions sum up together, Af, = + ^^iggsmo^ where the inter¬ 

esting one-loop contribution is given by [12-15] 


= -^l^MctanISlirh-f,^, (6a) 

^higgsino ^ -^^^tan 131 (6b) 
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There are also higher order calculations oi A}, available that are important for precision 
analyses [16-18]. 

The gluino loop contribution (6a) decouples with the gluino mass if the other SUSY 
parameters are fixed, but the higgsino one (6b) cannot be neglected for the desired values of 
jjL around the SUSY scale. For the numerical analysis in the course of this letter, we set Mg = 
Msusy which reduces for positive jJL. Moreover, we only include “active” third generation 
squarks as superpartners and implicitly take any other superpartner heavy (all gauginos 
besides the gluino which does not give a contribution to the effective Higgs potential at 
one-loop). 

There are handy exclusion limits, well-known for a long time, to simply check whether an 
unwanted, charge and color breaking (CCB) minimum appears for a given set of parameters 
in the MSSM. The constraints are on soft breaking trilinear couplings against soft breaking 
mass parameters as 


aI < 3(m22 + m^ -F m^). 


(7) 


see e. g. [19-26]. 

Mostly studied, however, are such couplings of up-type squarks to the up-type Higgs or of 
down-t3q3e sleptons to the down-type Higgs (where similar expression for down-type squarks 
can be obtained by relabeling the parameters). Couplings to the “wrong” Higgs doublet are 
mainly excluded in the analyses. The destabilizing contribution is always related to the 
trilinear part of the scalar potential, e. g. ~ ft has been shown [11] that the 

direction of the up-t3q3e Higgs field gets apparently destabilized from a (s)bottom loop effect. 
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In [11] only the field direction of the neutral Higgs, was considered—we now want to 
give a more complete view of the destabilizing effect leading to an analytic approximate 
exclusion on the combination ijlY^ in case the colored sbottom direction is included. Another 
exclusion can be obtained using a different direction in field space, where also the down-t 3 Ape 
Higgs scalar is needed. 

In this letter, we describe in the following section how to derive the analytic expression 
for the new CCB constraint from sbottom vevs and compare it to the numerical analysis of 
the global minima in the quantum (e. g. loop corrected) theory. Finally, we conclude. 

2 Finding CCB minima 

So far, we only discussed features of the scalar (one-loop) Higgs potential from Eq. (2) as 
described in [11]. In order to find the new (true) CCB vacuum, which hides behind the 
critical Higgs field value, we add to the potential of Eq. (2) (evaluated at = Mg^g^) the 
tree-level part of the sbottom potential, 

y-r = + Wfb, + bliml + \Y,hl\^)K 

- [blip-*Y^h°^ -A^hl)b^-\-h.c. ] + \Y^\^\b^\^\K\^ + D-terms. 

As was already pointed out before [27,28], the destabilizing term is always the trilinear 
one, so we expect a new stability condition for the combination juYj, taking 

Ab = 0. Actually, we cannot ignore D-terms in the tree-level potential accordingly to the 
neglect of all terms in the derivation of the one-loop Higgs potential, since also the 
Higgs self-couplings are ~ gi 2 - However, we can simplify (as usually done) the discussion 
considering so-called “D-fiat” directions. Those directions are most probably that kind of 
rays in field space in which unwanted minima develop. Non-D-flat directions are protected 
by the quartic terms that will always take over. The full D-term potential for the Higgs and 
sbottom scalar potential is given by 

2 ( 9 ) 

+- I'-SP+ 

We still ignore stop and stau fields and remark that the pure Higgs terms are already included 
in Eq. (2). Nevertheless, we make use of Eq. (9) to set the interesting directions; with b^ = 


5 


= b, we have the SU(3)(. D-flat direction. Considering the three-field scenario, we can 
reduce the degrees of freedom forcing all D-terms to vanish by the choice \h^\^ = |/i°p -I- |bp. 
Still rather large quartic terms survive in the potential, namely the |YbP terms from the F- 
term part in For that observation, we also look into a non-D-flat direction keeping 

^^ 8 ^^ + |bp)^, where the down-type Higgs is fixed at = v^j which is a constant and 

small number especially for large tan^,^ and therefore neglected with respect to potentially 
large field values of b and Note that contrary to most previous considerations [26,27,29, 
30] we are explicitly interested inb^O though A], = 0 and have tL,R = 0- In both ways we 
are considering a combined non-standard vacuum in the mixed sbottom and up-type Higgs 
direction instead of the pure down-down case. 

Let us figure out the analytic bound analogously to the famous A-parameter bounds like 
Uneq. (7), under which circumstances a CCB true vacuum appears. For that purpose, we 
shall choose the most probable field configuration that makes all the D-terms vanish. In the 
SU(3)(. X SU(2 )l X U(l)y D-flat direction, we assign b^^ = b^ = b and -I- b^. We 

consider only real fields and parameters now and in the following for simplicity. A different 
but not uninteresting bound will be derived in a direction where we keep the field strength 
at a fixed and small value, /i° = v^j 0. That way, we cannot reduce the quartic terms but 
still find a (new) analytic exclusion in the h° = b direction. 

= b An exact analytic derivation of the exclusion limits from the stability of the elec- 
troweak vacuum against formation of charge and color breaking minima is very easy to 
obtain in the one-field scenario. We follow the standard procedure which was pictorially 
reviewed in Ref. [8] . We collect the interesting parts of the tree-level potentials of Eqs. (2) 
and ( 8 ), 

= (M" - 2pLY^h)b^ + + A^b^ + + Xj,^h^b\ ( 10 ) 

with = tUq -\- m^, + pi^ and the self-couplings -F A,, = and 

X^b — This simplifies via b — h further to 

y= fh^h^ -Ah^-t Xh\ (11) 

^With V(j = V cos p we denote the standard electroweak vev of the down-type Higgs. 
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Figure 1: New exclusion limits including the formation of charge and color breaking minima with 
both (b) 7 ^ 0 and (b°) ^ 0. The left plot shows exclusions in the ju-tan^ plane with = —ISOOGeV 
which has been chosen conveniently to give the proper light Higgs mass within a few GeV and a com¬ 
mon soft breaking mass Mgusy = 1 TeV. The small points are excluded by the numerical comparison 
of the two minima (exclusion if CCB minimum deeper than trivial one). On the right-hand side we 
depict the crucial dependence of the non-standard minimum on the (mis) alignment of sbottom field 
and Higgs field value (b — ah with h = h^^, h^ = 0 and a e {0.8,1,1.1} for a given excluded point. 


with + m^, A = A;, + + Af,f, and A = We then find with the vev, 


V = (h) = 


3A+ y9A2 - 32m2A 
8 A 


and the requirement^ that for stable configurations V^jr, = > 0 , which is 

the new condition as (m^ is negative!) 

^u 




Y,^ + igl + gi)/2- 


( 12 ) 


Note that has a non-trivial dependence on tan^ and also A^ via Aj,, see Eqs. ( 6 ) 
and [12-15]. The (g^ + gl)/2 contribution is the left-over from the non-D-flatness which 
can be numerically of the same size as a threshold-resummed weakening the exclusion. 
This bound, however, does not fit exactly to the numerical exclusion as can be seen from 
Fig. 1 but provides an excellent approximation though actually (h) 7 ^ (b). The numerical 
exclusion limit shown in Fig. 1 agrees well with independent previous analyses on a similar 
situation [31] and are a bit stricter than the final results of [ 11 ], whereas a similar necessary 

^The potential of Eq. (11) reveals a strong first order phase transition, where the trivial minimum appears 
to be V(b = 0) = 0. Stable configurations need the potential value to be larger than that one. 
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condition was found for a slightly different direction in field space [32], 

|h°|2 = |h«|2 + |b|2 With the knowledge from above, it is straightforward to give a similar 
exclusion in the D-flat direction |/i°p = + \b\^. The remaining two-field scalar potential 

(real fields and parameters, = 0) can be further reduced aligning b — a/i° = ah with a 
(real) scaling parameter a: 

^D-flat = + ^22 ^ 2 m^ 2 '\/1 + h^ 

— 2yiY^a}b? + Y^(^2a^{l + a^) -I- 

that can be easily mapped on the expression of Eq. (11) resulting in the requirement that for 
stable configurations'^ 

2 2 

ml^(l + a^) + ml2±2ml2\/TT^+a^(m^^ + ml)> (14) 

This exclusion translated into the ju-tan^ plane is shown in Fig. 2 where we also display 
points that are excluded via the numerical minimization of the combined tree and one- 
loop effective potential. To enhance the significance of this bound (which is basically tan^- 
independent), we have employed running squark parameters in the tree-level sbottom poten¬ 
tial evaluated at the scale of the new minimum. Therefore, also corresponding parameters 
in the analytic exclusion (soft SUSY breaking masses and ju) have been taken at the same 
scale. Unfortunately, for the purpose of displaying the exclusion line, it is not clear at which 
scale those parameters have to be evaluated. As the second minimum genetically appears 
around one order of magnitude above the SUSY scale, we have set a fixed renormalization 
scale of lOMgusY therefore blue dots and the reddish area on the left-hand side Fig. 2 do 
not perfectly fit. Moreover, the excluded area by Uneq. (14) is not completely filled with ex¬ 
cluded blue points as there the sbottom-tree plus Higgs-one-loop potential shows a different 
behavior than the classical potential as also depicted in Fig. 3. 

Unequations like (14) or (12) follow from the tree-level potential and can be determined 
easily once a specific field line is selected. Going beyond tree-level changes the situation 
severely as can be seen from Fig. 3. A configuration which is obviously unstable (right-hand 
side) at the tree-level not even develops a second minimum considering the one-loop Higgs 
potential (the complete one-loop potential including sbottom directions was not employed 

'^The sign ambiguity origins from the fact, that we only need to constrain \h^\^ where the overall phase or 
sign is not constrained. 
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Figure 2: Exclusion in the /i-tan^ plane similar to the one shown in Fig. 1 (which is indicated by the 
grayish area) for the D-flat direction = |/t°p + \b\^. Blue dots have been excluded via numerical 
comparison of the two minima (if so) using the one-loop Higgs potential and an improved sbottom 
potential at the tree-level; the red line shows the exclusion of Uneq. (14) where the misalignment 
parameter a has been “fitted” for optical agreement of the blue dots and the reddish area to be 0.75; 
the actual a are different for each blue point. On the left-hand side, we have the —sign and on the 
right-hand side the -I—sign of Uneq. (14). 


for that purpose though should be available numerically). However, this effect is different 
in the “positive” direction where unstable configurations are driven towards more stable 
ones as can be seen from the left-hand side of Fig. 3. Usage of the renormalization group 
improved (tree-level) effective potential, where the couplings (Yukawa couplings and masses, 
actually no gauge couplings are they are absent in the genuine D-flat direction) are evaluated 
at a proper scale, ^ hint towards less restrictive exclusions. Where the tree-only potentials 
show non-trivial charge and color breaking minima, the loop-corrected potentials seem to 
stabilize the standard vacuum against formation of false vacua. 

Estimate of lifetime Are the developing charge and color breaking minima really a case 
for anxiety? As long as the lifetime of the “standard” electroweak vacuum is (much) longer 
than the present age of the universe, we basically do not have not worry and can take the 
issue of vacuum metastability for future generations. We estimate the lifetime of the desired 
vacuum for the scenarios provided in Figs. 1 and 3 using the triangle method of [34] and the 
instable potentials shown in the figures. However, similar to the scenario discussed in [11], 

^The choice of a proper renormalization scale is a bit vague and the decision whether to trust that choice 
in order to discard certain configurations is tenuous. For our purpose, we stick to the suggestion of Ref [33] 
and choose a scale Q = max as the largest field-dependent mass eigenvalue of the loop-contributing 

fields (in our case top and/or bottom (s)quark). 
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Figure 3: For a given parameter point (/i = 4TeV, tan^ = 40, Aj = l.STeV), we show exemplarily the 
behavior of the potential in the given direction in field space (and b = 0.75/t°). On the left-hand side, 

the positive sign for /t° was chosen, where the plot on the right has with real 

fields and parameters in both cases. The “tree + 1-loop” line means inclusion of the one-loop Higgs 
potential as of Eq. (2) plus the tree-level bottom squark potential (without D-terms each since they 
vanish by definition of the direction) evaluated with running parameters (soft-breaking squark masses 
and Yukawa couplings). For comparison, we show the “tree only” where the masses and couplings 
of the potential have been evaluated at the SUSY scale Mgusy — 1 TeV and the “RG-improved tree” 
potential where all soft masses and couplings are treated as running ones. 


where the decay time was found to be ridiculously small (details on the estimate have been 
given in [35]), we find our unstable solutions to be extremely short-lived concerning Fig. 1. 
This is not true for the genuine D-flat scenario shown in Fig. 3; here the lifetime is many 
orders the lifetime of the universe. 


3 Conclusions 

We have provided new (analytic) exclusion bounds in the MSSM from the formation of 
CCB minima. Contrary to previous considerations, we did not constrain the soft-breaking 
A-parameter by working in the direction of up or down fields only but connected the bot¬ 
tom squark direction with the up-type Higgs field. This procedure gives a constraint on 
where the bottom Yukawa coupling has an implicit dependence on the model parameters 
via Yj, = mi,/[Vd(l + A^')\. Under certain simplifications we have derived an analytic bound 
which is mostly in good agreement with the direct numerical exclusion from the minimiza¬ 
tion of the full (i. e. tree-level sbottom plus one-loop Higgs) effective potential considered 
in this letter. This bound complements existing CCB bounds and relates the bottom Yukawa 
coupling to soft SUSY breaking parameters (and the /r-parameter of the superpotential) 
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which is qualitatively different from existing traditional CCB bounds. The bottom Yukawa 
coupling itself depends nontrivially on the SUSY spectrum by virtue of threshold correc¬ 
tions for large tan^. A similar bound was found for the distinct direction in field space 
where all the D-terms vanish. The corresponding unstable solutions are rather metastable 
and very long-lived. Moreover, the comparison with quantum corrected potentials shows 
that even the metastable configurations tend to be stabilized by the loop contributions. This 
strengthens the previous bound in the explicit non-D-flat directions which stems from im¬ 
mensely short-lived configurations that persist in the presence of quantum corrections and 
is therefore more severe. The limitation to D-flat directions in the scalar potential as usually 
performed probably misses additional potentially dangerous directions. 

We constrained ourselves to cases with only one non-standard vev, accordingly the exclu¬ 
sions would change once more directions are taken into account. In those cases, however, 
the definition of flat directions suffers from ambiguities which makes the derivation of an 
analytic bound similar to Eq. (12) unclear. Similarly, the constraints can be extended to 
non-vanishing stop and stau vevs as has been done for the left-right mixing of staus [36] . 


Acknowledgments 

This work was supported by the Karlsruhe School of Elementary Particle and Astroparticle 
Physics: Science and Technology (KSETA). The author thanks U. Nierste for useful discussions 
on the topic and him and M. Spinrath for reading and commenting on the manuscript. 


References 

[1] S. R. Coleman and E. J. Weinberg, “Radiative Corrections as the Origin of Spontaneous 
S}Tnmetry Breaking”, P/ijs.Rev. D7 (1973) 1888-1910. 

[2] R. Jackiw, “Functional evaluation of the effective potential”, Phys.Rev. D9 (1974) 1686. 

[3] C. Ford, I. Jack, and D. R. T. Jones, “The Standard model effective potential at two loops”, 
Nucl. Phys. B387 (1992) 373-390, arXiv:hep-ph/0111190 [hep-ph] . [Erratum: Nucl. 
Phys.B504,551(1997)]. 

[4] S. R Martin, “Two loop effective potential for a general renormalizable theory and softly 
broken supersymmetry”, Phys. Rev. D65 (2002) 116003, 
arXiv:hep-ph/0111209 [hep-ph]. 

[5] S. E Martin, “Three-loop Standard Model effective potential at leading order in strong and top 
Yukawa couplings”, Phys. Rev. D89 no. 1, (2014) 013003, arXiv: 1310.7553 [hep-ph] . 


11 


[6] S. P Martin, “Four-loop Standard Model effective potential at leading order in QCD”, 

Phys. Rev. D92 no. 5, (2015) 054029, arXiv: 1508.00912 [hep-ph] . 

[7] J. Camargo-Molina, B. O’Leary, W. Porod, and F. Staub, “Vevacious: A Tool For Finding The 
Global Minima Of One-Loop Effective Potentials With Many Scalars”, 

Eur.Phys.J. C73 no. 10, (2013) 2588, arXiv: 1307.1477 [hep-ph] . 

[8] J. Camargo-Molina, B. O’Leary, W. Porod, and F. Staub, “Stability of the CMSSM against 
sfermion VEVs”, JHEP 1312 (2013) 103, arXiv: 1309.7212 [hep-ph] . 

[9] Y. Fujimoto, L. O’Raifeartaigh, and G. Parravicini, “Effective Potential for Nonconvex 
Potentials”, NucLPhys. B212 (1983) 268. 

[10] E. J. Weinberg and A.-q. Wu, “Understanding complex perturbative effective potentials”, 
Phys.Rev. D36 (1987) 2474. 

[11] M. Bobrowski, G. Chalons, W. G. Hollik, and U. Nierste, “Vacuum stability of the effective 
Higgs potential in the Minimal Supersymmetric Standard Model”, 

Phys.Rev. D90 no. 3, (2014) 035025, arXiv: 1407.2814 [hep-ph] . 

[12] L. J. Hall, R. Rattazzi, and U. Sarid, “The Top quark mass in supersymmetric SO(IO) 
unification”, Phys.Rev. D50 (1994) 7048-7065, arXiv:hep-ph/9306309 [hep-ph] . 

[13] M. S. Catena, M. Olechowski, S. Pokorski, and C. Wagner, “Electroweak symmetry breaking 
and bottom - top Yukawa unification”, NucLPhys. B426 (1994) 269-300, 

arXiv:hep-ph/9402253 [hep-ph] . 

[14] D. M. Pierce, J. A. Bagger, K. T. Matchev, and R.-j. Zhang, “Precision corrections in the 
minimal supersymmetric standard model”, NucLPhys. B491 (1997) 3-67, 

arXiv:hep-ph/9606211 [hep-ph]. 

[15] M. S. Catena, D. Garcia, U. Nierste, and C. E. Wagner, “Effective Lagrangian for the ibH~^ 
interaction in the MSSM and charged Higgs phenomenology”, 

NucLPhys. B577 (2000) 88-120, arXiv:hep-ph/9912516 [hep-ph] . 

[16] L. Hofer, U. Nierste, and D. Scherer, “Resummation of tan-beta-enhanced supersymmetric loop 
corrections beyond the decoupling limit”, JHEP 10 (2009) 081, 

arXiv:0907.5408 [hep-ph]. 

[17] A. Crivellin, L. Hofer, and J. Rosiek, “Complete resummation of chirally-enhanced loop-effects 
in the MSSM with non-minimal sources of flavor-violation”, JHEP 07 (2011) 017, 

arXiv:1103.4272 [hep-ph]. 

[18] A. Crivellin and C. Greub, “Two-loop supersymmetric QCD corrections to Higgs-quark-quark 
couplings in the generic MSSM”, Phys. Rev. D87 (2013) 015013, 

arXiv: 1210.7453 [hep-ph] . [Erratum: Phys. Rev.D87,079901(2013)]. 

[19] J. Frere, D. Jones, and S. Raby “Fermion Masses and Induction of the Weak Scale by 
Supergravity”, iVud.P/tys. B222 (1983) 11. 


12 


[20] M. Claudson, L. J. Hall, and I. Hinchliffe, “Low-Energy Supergravity: False Vacua and Vacuous 
Predictions”, NucLPhys. B228 (1983) 501. 

[21] J. R Derendinger and C. A. Savoy, “Quantum Effects and SU(2) x U(l) Breaking in 
Supergravity Gauge Theories”, Nucl. Phys. B237 (1984) 307. 

[22] C. Kounnas, A. Lahanas, D. V Nanopoulos, and M. Quiros, “Low-Energy Behavior of Realistic 
Locally Supersymmetric Grand Unified Theories”, NucLPhys. B236 (1984) 438. 

[23] L. Alvarez-Gaume, J. Polchinski, and M. B. Wise, “Minimal Low-Energy Supergravity”, 
NucLPhys. B221 (1983) 495. 

[24] L. E. Ibanez and G. Lopez, “N^l Supergravity, the Weak Scale and the Low-Energy Particle 
Spectrum”, NucLPhys. B233 (1984) 511. 

[25] M. Drees, M. Gluck, and K. Grassie, “A New Glass of False Vacua in Low-energy iV = 1 
Supergravity Theories”, Phys.Lett. B157 (1985) 164. 

[26] J. Casas, A. Lleyda, and C. Munoz, “Strong constraints on the parameter space of the MSSM 
from charge and color breaking minima”, NucLPhys. B471 (1996) 3-58, 

arXiv:hep-ph/9507294 [hep-ph] . 

[27] J. Gunion, H. Haber, and M. Sher, “Charge / Color Breaking Minima and a-Parameter Bounds 
in Supersymmetric Models”, NucLPhys. B306 (1988) 1. 

[28] R. Rattazzi and U. Sarid, “Large tan Beta in gauge mediated SUSY breaking models”, 

Nucl. Phys. B501 (1997) 297-331, arXiv: hep-ph/9612464 [hep-ph] . 

[29] A. Kusenko, R Langacker, and G. Segre, “Phase transitions and vacuum tunneling into charge 
and color breaking minima in the MSSM”, Phys.Rev. D54 (1996) 5824-5834, 

arXiv:hep-ph/9602414 [hep-ph]. 

[30] D. Chowdhury, R. M. Godbole, K. A. Mohan, and S. K. Vempati, “Charge and Color Breaking 
Constraints in MSSM after the Higgs Discovery at LHC”, JHEP 1402 (2014) 110, 

arXiv:1310.1932 [hep-ph]. 

[31] W. Altmannshofer, M. Catena, N. R. Shah, and E Yu, “Indirect Probes of the MSSM after the 
Higgs Discovery”, JHEP 01 (2013) 160, arXiv: 1211.1976 [hep-ph] . 

[32] W. Altmannshofer, C. Frugiuele, and R. Harnik, “Fermion Hierarchy from Sfermion Anarchy”, 
JHEP 1412 (2014) 180, arXiv:1409.2522 [hep-ph] . 

[33] G. Gamberini, G. Ridolfi, and F. Zwirner, “On Radiative Gauge Symmetry Breaking in the 
Minimal Supersymmetric Model”, NucLPhys. B331 (1990) 331-349. 

[34] M. J. Duncan and L. G. Jensen, “Exact tunneling solutions in scalar field theory”, 

Phys.Lett. B291 (1992) 109-114. 


13 


[35] W. G. Hollik, Neutrinos Meet Supersymmetry: Quantum Aspects of Neutrinophysics in 
Supersymmetric Theories. PhD thesis, Karlsruhe Institute of Technology, 2015. 
arXiv:1505.07764 [hep-ph] . 

[36] J. Hisano and S. Sugiyama, “Charge-breaking constraints on left-right mixing of stau’s”, 
Phys. Lett. B696 (2011) 92-96, arXiv: 1011.0260 [hep-ph] . [Erratum: Phys. 
Lett.B719,472(2013)]. 


14 


